Dipole emission and coherent transport in random media I 



M. DonairelB 

^ Departamento de Ftsica de la Materia Condensada, 
Universidad Autonoma de Madrid, E-28049 Madrid, Spain. 



o 
o 

(N 



O . 

"—I . 
> ' 

q 

O 

• I— I ' 

^ : 



> 
cn 

m 
o 



00 

o 

> 

X 
J3 



This is the first of a series of papers devoted to develop a microscopical approach to the dipole 
emission process and its relation to coherent transport in random media. In this Letter, we deduce 
general expressions for the decay rate of spontaneous emitters and the power emission of induced 
dipoles embedded in homogenous dielectric media. We derive formulae which apply generically to 
virtual cavities and, in the continuum approximation, to small real cavities. 

PACS numbers: 42.25.Dd,05.60.Cd,42.25.Bs,42.25.Hz,03.75.Nt 



It is a general issue in physics the characterization of 
a system by means of its coherent transport properties 
and the study of the decay of unstable local states. In the 
quite general case that the constituents of the medium 
couple to each other through dipole-dipole interactions, 
the decay of an excited particle takes place through radia- 
tive and non-radiative emission. In particular, it is known 
that the spontaneous emission rate, F, in a dielectric 
medium depends on the interaction of the emitter with 
the environment [ll]. This is so because the surrounding 
medium determines the number of channels into which 
the excited particle can decay. That is, the local density 
of states (LDOS). It is in this sense that the net effect 
of the random medium is to renormalize the vacuum as 
seen by the emitter. On the other hand, LDOS and F 
depend, upon additional properties of the emitter, on 
the statistical parameters which determine the coherent 
transport features of the medium. That is, on the elec- 
tric susceptibility x, the refraction index and the mean 
free path. As a matter of fact, null transmittance and 
inhibition of spontaneous emission are expected to occur 
in photonic band gap materials [3]. The understanding 
of life-times in random medium is relevant in the context 
of fluorescence biological imaging and nano-antennas 
On the other hand, understanding of unconventional 
coherent transport properties is essential in engineering 
metamaterials for electromagnetic and acoustic waves Q . 
In this Letter we deduce general formulae for the spon- 
taneous and induced dipole power emission in a homo- 
geneous random medium characterized by its electrical 
susceptibility. Longitudinal and transverse components 
are differentiated. It paves the way for the deduction of 
the relation between LDOS and F with coherent trans- 
port parameters. 

The power emitted in the process of spontaneous 
decay of a point dipole from an excited state ^{u) is di- 
rectly proportional to its decay rate F^. The relation is 
given by F^^ = ^W^, where 
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In the above equation, Tr is the trace operator, fj, is the 
dipole transition amplitude, lu is the transition frequency 
and (r , r) is the propagator of the field emitted by 
the dipole in the process back to itself. We denote 2- 
rank tensors with overlines. In general, /i, w and Guj all 
depend on the interaction between the emitter and the 
host medium. On the other hand, LDOS is proportional 
to 3|Tr{^<^(f,f)}| and so are F^; and [gj]. 

Next, consider that the emitter consists of the dipole 
induced by a fixed exciting field Eexdr) — -Eq (t^) on a 
scatterer of radius a and dielectric contrast ee(^)- ^ is 
assumed far from any resonance frequency of the emitter. 
Smallness implies a k^^ , kg = lu / c being the bare wave 
number. Thus, EQ{r) is uniform within the scatterer. If 
the emitter is one of the host scatterers in a homogeneous 
and isotropic random medium, the emitted power reads 
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where Xe — (^el^) — 1) -not to be confused with the 
susceptibility of the random medium- and 
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Qio) (^f'j being the tensor propagator of the electric field in 
free space. In spatial space representation, it consists of a 



Coulombian field propagator G^l (r) 



plus a radiation field propagator, G^,'^''^ (r) 
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In Eq.([2]), ^^{r,^) is the propaga- 
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tor of the field emitted at a point f' inside the particle 
back to another point f also within the particle. Notice 
that, because the emitter is polarizable in this case, the 
wave comes back and forth infinite times with propaga- 
tor Qi^. Hence, series Eq.®. In the above equation, 
for a <^ the electric field is nearly uniform and we 
can approximate / Q{r — a)Q^{r) ~ ^a^Guj{0)- This ex- 
pression is formally correct. However, the perturbative 
expansion of Guj (0) contains singularities hidden in both 
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the electrostatic and radiative parts of The later 

is related to resonance frequencies of the emitter polar- 
izability and so, negligible in the present case. The 
former is regularized by considering the finite size of the 
emitter, Lim{/ dV e(r - a)G^°'>ir)} = gi^I as a ^ 0. 
One way to avoid carrying this limit in further calcula- 
tions consists of dressing up the single particle suscepti- 
bility with all the 'in- vacuum' electrostatic corrections at 
once. That is, by defining Xe = j^f2^e ^^"^ electro- 
static polarizability ao = 47ra'^ T+2 • ^^"^ other hand, 
the part free of singularities can be written as 



■. =-t;(c_-l)0(r-a) 

m 

- = G'\ 
-»=G 



O = n 

(a) 



(*) 

2y^^ Tr{iliK<"''/j) ' 



-[l-i/3Tr{.0}r 



■■.[l-i/3Tr{0}]' 



FIG. 1: (a) Feynman's rules. (&) Diagrammatic representa- 
tion of Eq.Q. (c) Diagrammatic representation of the dress- 
ing up of Xe leading to ao. Approximation symbols denote 
that the field within the emitter is taken uniform, (d) Dia- 
grammatic representation of Eq.Q. 
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where are singularity- free and we have defined 

the II factors in the last row. In Eq.Q, the scripts || 
and _L denote, in Fourier space, the tensor components 
along and transverse to the propagation direction respec- 
tively. Note that longitudinal and transverse modes cou- 
ple to each other in the series Eq.© as scattering takes 
place off the dipole surface. With the above definitions 
we can rewrite Eq. ^ in terms of electrostatically renor- 
malized operators - compare with 0, Q , 
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where 
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The power emitted/absorbed by the induced dipole ac- 



cording to Eq.© is 
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The term in Eq.® corresponds to the power absorbed 
within the emitter. The term in Eq.© corresponds to 
the power radiated into the medium. The latter contains 
contributions of both coherent and incoherent radiation 
together with power absorbed by the host scatterers. 

Finally, consider the spontaneous emission of a point 
dipole like that in Eq.([T]), but now with an additional in- 
vacuum polarizability a = ao[l — ikQao/{6Tr)]~^. The sit- 
uation is analogous to that of a fluorescent particle placed 
on top of a host scatterer. Thus, the dipole moment of the 
system emitter-host-particle reads p = fl + ^aQ^{0)fl, 
a being the renormalized polarizbility of the host parti- 
cle. The spontaneous field emitted in the decay process 
gives rise to an induced dipole moment in the host par- 
ticle which modifies the decay rate of the emitter. It is 
plain from the perturbative development in Fig|2] that 
the emitted power can be written as 
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where we recognize again the power absorbed within 
the induced dipole in the last term and the power ra- 
diated into the surrounding medium in the remaining. 
Emission in vacuum is obtained by setting Gujir) — 

W ~ Im{(tI+i[G+0]-V«+«IOtQ]4>[O+i;;i}t»+...)i;.£] 

(a)' 

r'^~im{o([Q+ Q]4Q. [p. Q]-i,»[Q+ g^^m:-, ph..)-*^} 
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FIG. 2: (a) Diagrammatic representation of Eq.((2|. (6) Dia- 
grammatic representation of Egs. dQIlOp . 

^(0) (j.^ |;]2g equations above. 

For a complete description of W , F and LDOS we are 
just left with the computation of Gc,j{f,r) and its trace 
components, 27^ and 7|| in each case. Guj{r,r) is the 
propagator of the field throughout the bulk from the 
emitter back to itself. It includes multiple-scattering pro- 
cesses which are in general spatially correlated. Corre- 
lations are both those among host scatterers themselves 
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and those of host scatterers with the emitter. If the emit- 
ter does not perturb the statistical isotropy of the host 
medium, it occupies the center of a spherical cavity of 
radius R. On top of that, the infinite series in Eq.([3]) 
amount to multiple self-polarization processes. It is con- 
venient to formulate additional Feynman's rules to de- 
scribe the cavity and the self-polarization cycles. Those 
are given in FiglS] A self-polarization cycle carries a fac- 
tor [27_L-h7|i]5l = / d3re^^ff(r)(5(3)(f) in the perturbative 

expansions for W in place of {^)-^ J d^rg^''^r)e{r - 
a) -see FigUJ^c, d) . The cavity gives rise to a negative 
correlation function he — ~Q{r — R). 

In addition, we have to describe the field propagation 
throughout the bulk. Special attention is to be paid 
to the interaction between longitudinal and transverse 
modes. For simplicity let us assume the host medium is 
infinite. Therefore, no coupling to surface modes needs 
to be considered. A bulk propagator G"^ and dielectric 
and susceptibility tensors, and x'^ can be unambigu- 
ously defined. is the dyadic Green's hmction of the 
macroscopic Maxwell equations for the time-mode uj in 
the bulk. We drop the script oj hereafter in all quan- 
tities. X carries correlation effects and is thus made of 
one-particle-irreducible (IPI) multiple-scattering events. 
In Fourier space, translation invariance and isotropy al- 
low us to split Dyson equation for G{k) in two uncoupled 
and mutually orthogonal scalar algebraic equations, 



Gx(fc) = Gf (fc) - fco'Gf (fc)x±(fc)Gi(fc), (11) 
G||(fc) = g[|")(A) - fco^G|°)(fc)x||(fc)G||(fc), (12) 



where G^^^^k) 



(0), 



± vv = 4^ and Uf[k) = fe bemg 
a unitary vector along the propagation direction and 
A(fc) = I — k ® k being the projective tensor orthogo- 
nal to k. The longitudinal component is the propagator 
of the electrostatic field. Likewise, the transverse com- 
ponent is the propagator of the radiation field. In view 
of Eas. (|llll2p . longitudinal and transverse coherent pho- 
tons do not couple to each other as traveling throughout 
a random medium as it is the case of photons in free 
space. Eqs. (|llll2p can be solved independently yield- 
ing the rcnormalized propagator functions for the co- 
herent -macroscopic- electric field, G±{k) — — 



G||(fc) 



k®k 



_(k)-k^ ' 

However, a more detailed examen shows 



that longitudinal and transverse bare photons -i.e., nor- 
mal modes of G*^"^- do couple necessarily when they prop- 
agate in a random medium and experience multiple scat- 
tering processes. In Eas. (|llll2p . longitudinal and trans- 
verse bare photons enter both x±{k) and x\\ (k) by means 
of the spatial correlations among scatterers. In particu- 
lar, coupling between longitudinal and transverse modes 
shows up at the emitter cavity. We define the cavity 
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FIG. 3: (a) Feynman's rules. (6) Diagrammatic representa- 
tion of the bulk propagator G. 
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where cos 9 = k ■ k' . 

Consider first the emitter as either a spontaneous or 
induced dipole equivalent in all to the rest of host scat- 
terers or as a fluorescent particle on top of a host scat- 
terer. Therefore, it correlates to the surrounding as any 
host particle would do. The exclusion volume around 
the emitter is thus a virtual cavity as it does not per- 
turb the medium at all. In each self-polarization cy- 
cle, the emitted/polarizing field experiences in the bulk 
multiple-scattering processes which are correlated to the 
emitter /receiver. A typical process of these is depicted in 
FiglH^). There, only 2-point correlation functions, h{r), 
are used for simplicity. Note that, because the emitter 
and the receiver coincide, the correlations of any inter- 
mediate scattering process to one or the other extreme of 
the diagram are equivalent. This allows to attribute all 
the correlations to the emitter on the left, which amounts 
to the factor x/po; -FigslD^c, d). Therefore, we have 
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(2^)3 lpaki[l + x\\ik)] 
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where p is the density of host scatterers. The factor 2 
in front of 7_l stands for the two transverse polarizations 
while there is only one longitudinal. The term — 25R{7^''} 

accounts for the singularity of G^radi^)- This concludes 
the computation of the spontaneous emission of an emit- 
ter on top of a host scatterer as a function of the sus- 
ceptibility tensor x of the host medium and the scat- 
terer polarizability ao. We can split up the power emit- 
ted and the associated decay rates into longitudinal and 
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transverse components, 2rj_ and Fn, as it derives from 
Eq.® (compare with [^). It is plain by substitution 
of x(fc) ^ pa ~ {pa)'^klC in Eqs. (|llll2p that, for any 
hc(r) 7^ const, longitudinal and transverse photons do 
couple at the cavity surface. This is at the root of the con- 
tribution of longitudinal modes to the total decay rate, 
even in absence of absorbtion (see [13] )■ In particular, 
for i? <C fc(^^, the only contribution to Eqs. p3ll4|) comes 
from longitudinal photons -i.e. electrostatics, yielding 
C'li ±(fc) — This is recognized as the Lorentz-Lorenz 
(LL) cavity factor [10|. 

Next, let us consider that the emitter is either an 
impurity itself or seats on top of a polarizable impu- 
rity within an arbitrary cavity. That breaks manifestly 
translation invariance. Thus, the cavity is real and it 
is not strictly possible to define a Dyson propagator. 
The reason being that, in any multiple scattering pro- 
cess, beside the translation-invariant correlation func- 
tions joining host scatterers, each scatterer i is correlated 
to the emitter through the 2-point correlation function 
gc{fi) = 1 + hc{ri). A first approximation to go around 
this problem can be made if i? ^, ^ being the typical 
correlation length between host scatterers. In such a case 
we can consider that, for any IPI diagram in x(ri,r2), 
the emitter gets correlated to any of the host scatterers 
therein by simple convolution of x(n.7'^2) with gc{ri). 
This yields the series expansion depicted in Figini^6) in 
which wave propagation still depends on the emitter lo- 
cation. A further approximation consists of considering 
that, for fco^ ^ 1 the empty cavity do not affect the 
wave propagation but in its entrance and departure from 
the host medium. That is, both wave and emitter see a 
continuum beyond r ^ R. Thi s ap proximation is closer 
to the Onsager-Bottcher (OB) [13| approach as reformu- 
lated in Q (see also the quantum approach in Q). If 
G is the bulk propagator as defined in absence of cavity, 
the 7 functions take the form -see FigOJc), 
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+ KJ (|^[qi+G[|")]^G,|xii(fc), 



(18) 



where an analogous equivalence to that illustrated in 
Fig|3{6) has been used. 

In summary, we have derived general expressions 
for the longitudinal and transverse components of the 
power emission and the decay rate of induced dipoles 
and spontaneous emitters respectively in function of 
the electrical susceptibility of the host medium. These 
are Eas. (|ll7l8l9ll0p . depending on the emitter nature, 
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FIG. 4: (a) 2-point correlation function h(r) and associated 
interaction vertex, (b) Diagrammatic representation of the 
equivalence between multiple-scattering processes amounting 
to Q. (c),(d) Alternative representation to Fig[ljd). 
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A ± /I ft' ir^ij, ft i I V 



FIG. 5: (a) 2-point correlation function gdj). (b) First ap- 
proximation of W^. (c) Diagrammatic representation of the 
approximate formulas in Eqs. p7ll8|l . 



together with the appropriate 7-factors. Those are, 
Eas. (ll5ll6p for a source seated in a virtual cavity. They 
arc exact and their diagrammatic representation is that 
of Figl2j6) -up to constant prefactors. For a source 
seated in a real cavity drilled in a continuous medium, 
the 7-factors are in Eqs. p7ll8p . They are approximate 
and their diagrammatic representation is that of Fig[5l[c). 
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